Astronomy & Astrophysics manuscript no. Artigol 'DC"structabstract2 


©ESO 2012 


December 13, 2012 





Impact of stellar companions on precise radial velocities 

D. Cunha''2, p. Figueira', N. C. Santos' ^ c. Lovis^ and G. Boue'-'' 



O 
o 

o 

<N 



Oh 

Oh! 

o 



> 

00 

00 
(N 

(N 
(N 



X 



Centro de Astrofisica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 
e-mail: Diana . Cunha@astro . up . pt 

Departamento de Fisica e Astronomia, Faculdade de Ciencias, Universidade do Porto, Rua do Campo Alegre 687, 4169-007 Porto, 
Portugal 

Observatoire Astronomique de I'Universite de Geneve, 51 Ch. des Maillettes, - Sauvemy - CH1290, Versoix, Suisse 
ASD, IMCCE-CNRS UMR8028, Observatoire de Paris, UPMC, 77 avenue Denfert-Rochereau, 75014 Paris, France 



Received December 13, 2012; accepted November 29, 2012 

ABSTRACT 

Context. With the announced arrival of instruments such as ESPRESSO one can expect that several systematic noise sources on the 
measurement of precise radial velocity will become the limiting factor instead of photon noise. A stellar companion within the fiber 
is such a possible noise source. 

Aims. With this work we aim at characterizing the impact of a stellar companion within the fiber to radial velocity measurements 
made by fiber-fed spectrographs. We consider the contaminant star either to be part of a binary system whose primary star is the target 
star, or as a background/foreground star 

Methods. To carry out our study, we used HARPS spectra, co-added the target with contaminant spectra, and then compared the 
resulting radial velocity with that obtained from the original target spectrum. We repeated this procedure and used diff'erent tunable 
knobs to reproduce the previously mentioned scenarios. 

Results. We find that the impact on the radial velocity calculation is a function of the diff'erence between individual radial velocities, 
of the diff'erence between target and contaminant magnitude, and also of their spectral types. For the worst-case scenario in which 
both target and contaminant star are well centered on the fiber, the maximum contamination for a G or K star may be higher than 10 
cm/s, on average, if the diff'erence between target and contaminant magnitude is Am < 10, and higher than 1 m/s if Am < 8. If the 
target star is of spectral type M, Am < 8 produces the same contamination of 10 cm/s, and a contamination may be higher than 1 m/s 
if Am < 6. 

Key words. Planets and satellites: detection - Techniques: radial velocities 



1. Introduction 

The search for extrasolar planets is currently a very active 
field of research in as tronomy. Since the first discovery of 
iMayor & Ouelozl (Il995l) . more than 800 planets were discov- 
eredQ, of which 80 % were detected with the radial velocity 
method (RV). Because it is the workhorse for planetary detec- 
tions, this method received much attention from the community, 
which const antly increased its p recision (for some of the latest 
results check lMayor et al.ll201 ll) and tried to characterize its lim- 
itations. 

One of the fundamental drawbacks is that it is an indirect 
method, and therefore one should be extremely careful with false 
positives, i.e., signals created by other planetary RV signals. The 
literature provides many examples and techniques used to pin- 
point these signals (e.g. iQueloz et a l."200l': 'Santos et al"]l2002l: 
iMelo et al.|[2007l: iHuelamo et al.l 120 08: Fisueira et aLl|2010|) or 
to average them out (iDumusque et al.i i201 Ibua) . With the ad- 
vent of very high precision s pectrographs, such as ESPRESSO 
or CODEX (iPepe et al.ll2010l) , the RV precision enters a new do- 
main, that of cm/s, which in turn will increase the quality of our 
characterization o f these unwanted yet ubiquitous signals (e.g. 
ICegla et al.l[20T2h . 

Here we address a so far unexplored mechanism that is capa- 
ble of distorting RV signals and creating false ones: the contam- 
ination of the main spectrum by that of a (usually unresolved or 

1 ,http://exoplanet.eu/ as of 27/09/2012 



otherwise undetectable) companion. We consider two diflTerent 
cases: 

1 . that of a faint gravitationally bounded companion and 

2. that of an unbound star that is aligned with the target at the 
moment of observation. 

We assumed the target stars to be placed at the Galaxy disk edge 
to consider a representative case that encompasses several sce- 
narios of contaminant stars. We used real high signal to noise ra- 
tio (S/N) observations obtained with HARPS to create the com- 
posite spectra and process these with the RV calculation pipeline 
to evaluate the impact of the contamination in the most realistic 
way possible. While we aim at characterizing this eflTect down to 
the precision level of ESPRESSO, 10 cm/s, it can be assessed at 
an even lower precision level. 

Although this study is performed to characterize the impact 
on the RV when using the cross-correlation function (CCE), 
point spread function (PSE) modeling will also be aflTected by 
stellar companions, if in a diflTerent way. The results of this study 
are particularly interesting for deep surveys in crowded fields 
(Kepler). 

If the contaminant RV remains constant, stellar blends will 
only cause a undetectable contamination in the RV calculation. 
However, if the contaminant RV is changing, if the seeing varies, 
or if there are guiding problems during observations, the eflfect of 
the contamination will no longer be constant and it might mimic 
the presence of a planet. As an example of a false positive caused 
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by blended stars we refer to the work of ISantos et al.l (l2002l) . in 
which RV measurements derived from CORALIE blended spec- 
tra of HD 41004 AB have unveiled a radial- velocity variation 
with a period of 1.3 days and a small amplitude of 50 m/s, com- 
patible with the signal expected to be caused by a planetary com- 
panion to HD 41004 A. Another example is WASP-9b, which 
was discarded as a planet after it was discovered that the signal 
was due to a fortuitous alignment (Triaud, priv. comm.). 

In Sect. [21 we explain which stars we consider most likely 
to contribute to a shift in the RV calculation, if the star is ei- 
ther in a G-, K- or M-dwarf binary system, or is affected by a 
fortuitous alignment. Then, in Sect. [3] we describe the properties 
of the spectra used in our simulations. In Sect. |4] we describe 
the method used to scramble spectra from target star and stel- 
lar companion. Results from our work, assuming the worst-case 
scenario in which both target and contaminant are well centered 
on the fiber, are shown in Sect. [5] The results are followed by a 
statistical analysis in Sect. 15.41 which is based on these results, 
but assumes a more realistic situation, in which the contaminant 
star may not be fiber centered. In Sect. [6] we discuss our work, 
and we present our conclusions in Sect. [7] 

2. Contaminant cases 

When one observes a target star, a secondary star may be present, 
even without one's knowledge. This secondary star may be a 
source of contamination for the calculation of the target star's 
RV. Secondary/contaminants stars may appear as two diff'erent 
types: Real binaries or fortuitous alignments. 

2.1. Real binaries 

In the case of a real binary system, the contaminant star will be 
gravitationally bound. In this section we present the probability 
distributions of real binary properties that we used in our work. 

He re, for G- and K primaries we based our study on the re- 
sults of iDuquennoy & Mayod (Il99ll) . who studied the multiplic- 
ity among solar- type stars in the solar neighborhood. They con- 
sidered a subsample of 164 primary stars with spectral types be- 
tween F7 and G9. Using Table 7 of Duqu ennoy & Mayo r ( 199i|) 
on the mass-ratio distribution among G-dwarfs binaries, we can 
calculate the probability of observing a binary with a mass-ratio 
q = M2/M1, where M2 is the mass of the secondary star and Mi 
the mass of the primary. These probabilities are shown in Table 
[T] in which we can see that there is a maximum for the mass- 
ratio range ]0.2,0.3]. The distribution of orbital periods (P) of 
their study is approximated to a Gaussian with logP = 4.8 and 
cTiogP = 2.3. 

Duquennoy & Mayor also presented the eccentricity distri- 
bution for P < 1000 and P > 1000 days, shown here in Table 

m 

We stress that there is a large uncertainty for VLMC (very 
low mass companion) binaries (with ^ < 0.1), and that this study 
does not include values of q below 0.01. We also emphasize that 
this mass ratio distribution is calculated for binary stars, but we 
recall that, also according to Duquennoy & Mayo3 (ll99lb . there 
are 57% of binary systems with a mass ratio higher than 0.1, and 
43% of apparently single stars. Of these 43%, (8 ± 6)% most 
probably have a VLMC, which puts the percentage of single 
stars at ~ 30%. 

For M-dwarfs we used the distrib utions of mass ratio , q, and 
semi-major axis, a, from the work of iJanson et al. I (I20T1 on M- 
dwarf multiplicity. The authors suggest an M-dwarf total multi- 



Table 1 : Probability of observing a binary with a mass ratio q = 
M2/M1 amon g G-dwarf binaries for perio ds P < 10"^ days and 
P < 10^ (from lDuquennov & Mavo3ll99lh . 



^max 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 


P(q) [%] 


P < Wdays 

12.0 16.3 13.1 12.7 12.3 12.9 6.5 1.9 2.6 3.4 6.2 


P(q) [%] 


P > lO^'days 

11.1 11.9 18.4 15.8 10.5 11.2 6.6 6.6 4.0 4.0 


P(q) [%] 


Total 

11.4 13.6 16.4 14.6 11.2 11.8 6.5 4.8 3.4 3.8 2.4 



Table 2: Probability of observing a binary with an eccentricity e 
among G-dwarfs binaries for periods P < 10^ days and P < 10^ 
(from iDuquennoy & Mayor. 1 99 li) . 



^max 


0.15 


0.3 


0.45 0.6 0.75 


0.9 1 


P(e) [%] 


12.50 


43.75 


P < lO^days 
31.25 6.25 6.25 


0.0 0.0 


P(e) [%] 


5.88 


11.76 


P > lO^days 
23.53 14.71 20.59 


23.53 0.0 



Table 3: Probability of observing an M-dw arf binary with a mass 
ratio q = M2/M1 (from lJanson et al.ll2012h . 



^max 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


P(q) [%] 


0.0 


1.61 


5.95 


10.45 


12.93 


13.84 


14.48 


13.03 


14.58 


13.13 



Table 4: Probab ility of observing a n M-dwarf binary with a cer- 
tain log a (from lJanson et al.|[2012h . 





0.25 


0.5 


0.75 1. 


1.25 


1.5 


1.75 


2. 2.25 


2.5 


Pdogfl) [%] 


0.51 


2.60 


9.02 15.39 


17.13 


17.65 


14.73 


11.74 7.85 


3.37 



plicity fraction of 34% . Table[3]shows the probability of observ- 
ing an M-dwarf binary with a mass ratio q. Janson et al. found 
a uniform distribution to be more consistent with the mass ratio 
of their observed sample, than a rising distribution. Accordingly, 
we assumed their unifor m distribution in our work. For the semi- 
major axis distribution iJanson et"an (1201 2l) also compared the 
Sun-like and the narrow distribution, finding the latter to be more 
consistent with their sample. We therefore used this distribution, 
which is shown in Table IH As no information was given on the 
eccentricity distr ibution for M-dwarf b i naries , we assumed the 
distribution from IDuquennoy & Mayoj (Il99lb . 

These distributions were used in our statistical analysis 
(Sect. El. 

2.2. Fortuitous alignment 

Even if the target star is a single star, an alignment with a back- 
ground/foreground object may occur. To study the probability 
of sii ch an event we used the Besangon model (iRobin et al.l 
I2003h . This allows one to compute the probable stellar con- 
tent on a given direction of the Galaxy, permitting one to in- 
fer the most probable contaminants. To do so, we simulated 
an observation in the direction of the HD 85512 coordinates. 
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Table 5: Most probable contaminants for the fortuitous align- 
ments. 



Spectral type 




P [%] 


[^0, F51 


[13, 16] A [20,23] 


0.08 


[F5, G0[ 


[13,24] 


1.25 


[GO, G5[ 


[13,25] 


1.74 


[G5, i^0[ 


[13,25] 


1.32 


[KO, K5[ 


[12,26] 


2.5 


[K5, M0[ 


[16,30] 


4.48 


[MO, M5[ 


[17,30] 


16.99 


[M5,M8[ 


[20, 30] 


12.12 


[M8, M9] 


[21,30] 


13.17 



Table 6: Stars used in our simulations, their spectral type, visual 
magnitude (my) and S/N at the center of the spectral order 50 
(SN50). 



Object 


Spec.Type 




SN50 


HD20852 


F2 


6.92 


225; 98 


HD103774 


F5 


7.12 


102; 146 


Sun 


G2 


3.65 


225; 260 


Tau Ceti 


G8 


3.50 


300; 450 


HD69830 


KO 


5.95 


253; 380 


Alpha Cen B 


Kl 


1.33 


385; 399; 302 


HD85512 


K5 


7.65 


156; 86 


G1436 


Ml 


10.68 


15; 31; 46 


G1581 


M3 


10.57 


37; 22 



from Simbad (http:simbad.u-strasbg.fr/simbad/sim-fid) 
rriv of the moon within the HARPS 0.5" radius fiber. 

(/, b) = (271.6759, 8.1599), which is positioned at the border of 
the Galactic disk. We also considered a distance interval of 50 
kpc (approximately the diameter of the IMilky Way), magnitudes 
between and 30, and a solid angle for a radius of 3°. The objec- 
tive was to have a large number of stars for each bin of magni- 
tude and spectral type and then to divide it by the ratio between 
the area considered and that for the HARPS fiber (radius of 0.5 
arcsec), thus obtaining the density of stars or the probability to 
have a fortuitous alignment within the HARPS fiber. In Fig.[Tlwe 
present a contour plot of the total number of stars obtained in the 
simulation for each bin of magnitude and spectral type. For more 
details on these values and the values for the star density when 
normalized for a HARPS size fiber, see Appendix [A] The most 
probable contaminants for the accidental alignments are shown 
in Table El 

3. Data 

We contemplated target stars from spectral type G to late M, 
with visual magnitudes, m^, between 5 and 15. As described in 
the previous section (Sect. 12.21) , the most probable contaminant 
stars are of spectral types FGKM. To represent both target and 
contaminant spectral types, we chose spectra from nine stars, 
which are shown in Tabled 

We used spectra with an S/N at the center of the spectral 
order number 50 (which corresponds to a wavelength of 437.28 
nm), varying between 15 for an M star, and 450 for a G star. The 
S/N values of these spectra can be found in Table [6l The stars 
were chosen because they had the highest S/N in their spectral 
type range, as observed with HARPS. 



4. Method 

The goal of this work is to study possible stellar contamination 
of a target star. To do so, we co-added target and contaminant 
spectra. We not only considered contaminant stars of spectral 
types F2, F5, G2, G8, KO/Kl, K5, Ml, and M3, but also diff'erent 
ratios of magnitude between target star and contaminant as well 
as diflTerences in RV. To sum two spectra with these particular 
conditions is not trivial. We need to control the magnitudes and 
RVs. 

We set the magnitude V of the target (ntt) star according to 
the data available in SimbadH and then calculated the magnitude 
of the contaminant (m^) with respect to the target star using the 
apparent magnitude - flux relation. We also considered the rel- 
ative RV between the two spectra as a tunable knob. However, 
when processing our spectra using the HARPS pipeline the re- 
duction takes into account the header (and properties) of the par- 
ent spectra. To overcome this situation we calculated the contam- 
inant RV in the target reference frame and shifted it afterwards. 

We then co-added the target spectrum with the modified con- 
taminant spectra and processed the result with the HARPS RV 
pipeline. 

More details of the method are presented in Appendix |B] 

4.1. The "real binaries" star strategy 

To study the impact of a companion on the RV calculation of 
a binary primary star, we reproduced the cases discussed in the 
work of iDuquennoy & Mayoj (Il99ll) . We considered a primary 
G2 star with six diflTerent secondary/contaminant stars: G2, G8, 
Kl,K5,Ml,and M3. 

The first case had a G2 as contaminant with positive and neg- 
ative values of ARV = RVC^ - RVCu where RVCt is the target 
bary centric radial velocity (drift corrected) without contamina- 
tion, and RVC^ is the contaminant barycentric radial velocity 
(drift corrected) in the target reference frame. In Fig. 0we show 
the impact, RVCt - RVC'^, of the stellar companion on the target 
RV, where RVC^ is the target barycentric radial velocity (drift 
corrected) with contamination. 

For each target-contaminant pair we considered a range of 
\ARV\ between and 30[Km/s] (in steps of 1 [Km/s]), except 
for a secondary star of spectral type K5 in which the maximum 
\ARV\ was of 24 [Km/s] (see Appendix |B] for details). 

4.2. Fortuitous alignment strategy 

The fortuitous alignment case is more generic. Here we consid- 
ered six spectral types of primary stars (G2, G8, KO, K5, Ml, 
and M3) in combination with the eight spectral types as contam- 
inants - the same six used as target stars plus an F2 (HD20852) 
and an F5 (HD 103774) star. For each spectral type of primary 
stars we considered contaminants with magnitudes mv between 
5 and 15 (in steps of one). The contaminants considered were 
those described in Table [5] 

Although the Besangon model includes oxygen and carbon 
AGB (asymptotic giant branch) stars and white dwarfs in the 
search, we only considered main-sequence (MS) stars. Very 
metal-poor stars were also excluded from our study. For an Ml 
target star a KO star (HD69680) was used instead of the Kl men- 
tioned above. This was because the barycentric Earth radial ve- 
locity (BERV) value form the Kl star did not allow the desired 
shift in RV For each target-contaminant pair we considered the 

^ http://simbad.u-strasbg.fr/simbad/sim-fid 
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[o8,05[ [O5,B0[ [B0,B5[ [B5,A0[ [A0,A5[ [A5,F0[ [F0,F5[ [F5,G0[ [G0,G5[ [G5,K0[ [K0,K5[ [K5,M0[ [M0,M5[ [M5,M8[ [M8,M9] 

Spectral type 

Fig. 1: Number of stars in each bin of spectral type and magnitude (radius of 3°). 

5. Results 

5.1. Real binaries 

The results for the contamination by the secondary star of a bi- 
nary with a primary star of spectral type G2 as a function of ARV 
can be found in Fig.O Obviously, the impact on the RV calcula- 
tion is a function of the difference between target and contami- 
nant RV (ARV). One can interpret a CCF as an "average" line of 
a spectrum, that is fitted by a Gaussian curve to measure the cen- 
ter of the CCF. In the presence of a second star, the CCF is com- 
posed of two curves. For ARV = the two curves fully overlap, 
and so the impact on the Gaussian fit and on the RV calculation 
is minimum. As ARV increases, the impact on the Gaussian fit 
increases until the point when the two CCF curves overlap, with 
the contaminant RV, \RVC[.l close to 0.6 times the FWHM of the 
target CCF. When the CCF curves become distinct, the impact on 
the measured RV of the primary decreases. For \ARV\ > 20km/s 
the peak of the contaminant CCF curve exits the CCF window 
that is usually computed with a value of 20 km/s, and for values 
\ARV\ > 25 km I s the contaminant CCF curve completely exits 
the target CCF window. Nevertheless, the contamination is not 
null. There are bumps in the wings that slightly contaminate the 
target spectra. 

Figure |4] shows the maximum contamination for each type 
of contaminant. As we can see, as the secondary stars become 
fainter and of later spectral type, the maximum contamination 
decreases. 

5.2. Contamination by the sky 

The sky brightness is not zero and its light results from the re- 
flection and scattering of Sun spectra. The sky brightness has a 



Target G2 (mv=3.65) with contaminant G2(m\/=3.65) 



-20 -10 10 20 

|RVCc'-RVCt| (km/s) 

Fig. 2: Shift of a G2 target RV caused by a contaminant of the 
same spectral type and magnitude, my, in function of the difl'er- 
ence between individual radial velocities. 



difl'erences between its RV, |A7?y|, under 20 Km/s, which is half 
the width of the CCF window. Although giant stars will often be 
contaminant stars, we assumed that spectra from MS stars are 
similar to those of giants and used only spectra of MS stars in 
our simulations. 
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contaminant G2 




12 14 16 18 20 22 24 



contaminant Kl 




12 14 16 18 20 22 24 



m m f • • 



contaminant Ml 

-| \ 1 r 




12 14 16 18 20 22 24 



• • - ■ • • • 
• • • • ■ 



2000 
1500 
1000 
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80 
60 
40 
20 

3.0 

2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



10 15 20 25 

|RVCc'-RVCf| (km/s) 



contaminant G8 




12 14 16 18 20 22 24 



contaminant K5 




12 14 16 18 20 22 24 



contaminant M3 

-| r 




12 14 16 18 20 22 



10 



15 



20 



25 



Fig. 3: Impact on the RV of a G2 target star caused by a companion of spectral type G2, G8, Kl, K5, Ml, and M3 (left to right, top 
to bottom). For each case a close-up with a bump for the impact for 10 < ARV < 24 km/s is shown. 



magnitude per arcsec^ in V band of 21.8 on a new-moon night, 
and of 20.0 on a full-moon night. To study the impact of this con- 
tamination we assumed that sky spectrum is the solar one. This 
is equivalent to say that for a fiber of 0.5" radius, we have a sys- 
tematic contamination of a G2 star with my ~ 20 on a full-moon 
night to a limit of a G2 star with my ~ 22 (new moon). 

The results of the maximum contamination by the sky bright- 
ness study are shown in Fig. [5] We can see that when observing 
a G2, G8, KG, or a K5 target star fainter that my = 11, or Ml, 
M3 fainter that my =14, the errors induced on the observed RV 
can be larger than 10 cm/s, which can be a problem for future 
instruments such as ESPRESSO. For a night of observation with 
a full moon, induced errors become larger than 10 cm/s for G- 
and K target stars fainter that my = 9 and for an M star fainter 
than my =12. 

5.3. Fortuitous alignment 

In this section we present a more exhaustive study of all pos- 
sible sources of stellar contamination. Beyond the moon as a 
contaminator, there is a vast range of possible contaminators, as 
was discussed in Sect. 12.21 Our study is simplified because our 
analysis holds for similar spectra with diff'erent magnitude for 
target and contaminant star, but with same diff'erence between 



them (Am = - m^), as we can see in Fig.O This way, we only 
carried out our study on the impact for a target star in which my 
diff'ers \Am\ from the contaminant star. The result of our study 
for contamination as a function of spectral type is shown in Fig. 

m 

We first note that there is a clearly defined relation between 
Am and the maximum of the contamination, except for some 
cases where the diff'erence in magnitude between target and con- 
taminant is null. A misidentification of the companion correla- 
tion peak with that of the primary may happen due to a (slightly) 
larger contrast. If we take only Am > 1, we can describe the 
relation between Am and the maximum of the contamination as 

M2ix\RVCt - RVQ = 10^^'^+^ (1) 

where a and b are the coefficients presented in Table |71 In this 
table we also present the value for the root mean square (rms) 
for each fit. By comparing the panels of Fig. [7] we can also 
see that the contamination depends not only on Am, but also 
on the target-contaminant spectral type combination. This be- 
comes more evident in Fig. [8l This figure shows the maximum 
contamination as a function of the contaminant spectral type for 
a G2 and M3 target star with Am = 10. We find two reasons 
for this behavior. First, less contamination occurs with increas- 
ingly different spectra. Also, the contamination will depend on 
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Fig. 5: Maximum contamination of the sky brightness on a a) G2, b) G8, c) KO, d) K5, e) Ml, and f) M3 target star spectra. The 
blue triangles show the impact on the RV caused by the sky brightness on a full-moon night, while black diamonds show the impact 
on a new-moon night. 



the depth of spectral line. Thus, for similar spectra, cooler stars 
will cause a higher contamination, because their spectral lines 
depth is higher. 



5.4. Statistical analysis of the contamination 

After calculating the impact caused by a contaminant star in the 
RV calculation, it is important to known how often this occurs. 
Therefore we calculated the probability of measuring the target 
RV with a given contamination. 

In the previous two subsections (Sects. [5Jl and 15.31) we pre- 
sented a study for the impact to have a contaminant star within 
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Fig. 6: Contour plot of the maximum induced RV shift in the tar- 
get \RVCt-RVCf \ depending on the target (m^) and contaminant 
(mc) magnitudes. 



the fiber on the RV calculations, assuming that the fiber was 
well centered on target and contaminant. Although this is a rel- 
evant study, the most probable scenario for real observations is 
that of a contaminant star that not centered in the fiber. In this 
case, we cannot consider that both target and contaminant fluxes 
will completely enter the fiber, and therefore we calculated how 
much of each flux enters the fiber. The fraction of binary and for- 
tuitous alignments and their properties also follows established 
probability laws. To take these parameters into account in a real- 
istic way, we ran Monte Carlo simulations with 10000 trials, to 
calculate the distribution function of the impact of the presence 
of a contaminant star on the RV calculation. Once again, we had 
to consider both cases of contamination: binaries and fortuitous 
alignment. If in our simulations one star was in a binary sys- 
tem and had a background/foreground at the same time, we only 
considered the stronger RV eff'ect. In these simulations we did 
not take into account the contamination by the sky brightness 



Table 7: Coefficients a slxi b for each fitting of the relation be- 
tween Am and the maximum impact on the RV calculations: 
Msix\RVCt - RVCf\ = lO^-^'^^^, and rms resulting from each 
combination of target spctral type (TST) and contaminant spec- 
tral type (CST). 



TST 


CST 


a 


b 


rms 
[cm/s] 


G2 


F2 


-0.4008 


5.0733 


815.0813 




F5 


-0.4019 


5.1601 


1644.9045 




G2 


-0.4022 


5.4288 


4491.1375 




G8 


-0.4023 


5.4829 


4956.4768 




Kl 


-0.4026 


5.4556 


5136.0605 




K5 


-0.4009 


5.1859 


1500.0546 




Ml 


-0.4001 


4.8096 


23.3793 




M3 


-0.3996 


4.7301 


197.1604 


G8 


F2 


-0.4008 


4.9379 


681.0614 




F5 


-0.4019 


5.0222 


1075.0923 




G2 


-0.4027 


5.3440 


3968.9098 




G8 


-0.4021 


5.3154 


2697.4956 




Kl 


-0.4027 


5.3641 


4215.8544 




K5 


-0.4007 


5.11 


1139.4446 




Ml 


-0.3995 


4.7357 


3.0972 




M3 


-0.3984 


4.651 


51.5545 


KO 


F2 


-0.3812 


4.7320 


1566.6309 




F5 


-0.4008 


4.9521 


871.3506 




G2 


-0.3933 


5.2712 


5307.563 




G8 


-0.3929 


5.2613 


4246.1257 




Kl 


-0.3954 


5.3704 


7390.8788 




K5 


-0.3914 


5.1827 


3347.1766 




Ml 


-0.3799 


4.7540 


566.5518 




M3 


-0.3777 


4.6783 


532.1239 


K5 


F2 


-0.3988 


4.9559 


200.036 




F5 


-0.3999 


5.0865 


724.4446 




G2 


-0.4035 


5.4592 


7366.4185 




G8 


-0.4032 


5.5354 


8587.4134 




Kl 


-0.4053 


5.5580 


14146.5377 




K5 


-0.4030 


5.3793 


4935.42 




Ml 


-0.4011 


4.9458 


71.0313 




M3 


-0.4009 


4.9257 


98.6394 


Ml 


F2 


-0.4056 


3.5177 


18.1497 




F5 


-0.3932 


3.4673 


24.6295 




G2 


-0.4018 


4.1212 


10.1997 




G8 


-0.4011 


4.2786 


13.0058 




Kl 


-0.3990 


4.1604 


24.4166 




K5 


-0.4014 


4.5525 


19.9346 




Ml 


-0.3985 


5.0993 


33.3093 




M3 


-0.3994 


5.1581 


56.2493 


M3 


F2 


-0.4052 


3.4770 


14.7829 




F5 


-0.3924 


3.4421 


25.7393 




G2 


-0.4035 


4.0854 


30.8372 




G8 


-0.4019 


4.2469 


23.5459 




Kl 


-0.4024 


4.2649 


35.7392 




K5 


-0.4021 


4.5317 


32.2298 




Ml 


-0.3991 


5.0440 


13.2551 




M3 


-0.3999 


5.1436 


63.9946 



either, and also just considered that only ARV < 20 Km /s will 
produce a significant impact. We ran the simulations for each 
case of spectral type target considered in previous sections, at a 
distance d = 1, 5, 10, 50, and 100 [pc] from Earth. 

Before presenting our results of the statistical analysis we 
describe in more detail t he con siderations to statistically analyze 
the real binaries (Sect. 15.4.11 ) and fortuitous alignments (Sect. 
[5A2|. 
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Fig. 8: Maximum contamination \RVCt - RVCf \ as a function of 
the contaminant spectral type for a G2 (blue dots) and M3 (green 
diamonds) target star with Am =10. 



5.4.1 . Statistical analysis of binaries 

For our statistical analysis of contamination in binary systems, 
we considered the mass ratio distributions , q, mentioned in Sect. 
[2T|_We used Table[T]based on the study of iDuquennoy & Mayod 
(Il99ll) for target sta rs of spectral types G and K, and Table [3] 
based on the study of lJanson et al.l (l2012h for target stars of spec- 
tral type M. Because we had Ml and M3 spectra for M stars, we 
used the Ml spectra for [MO, M4[ stars, an d for the later types we 
used t he M3 spectrum. From the work of IDuquennoy & Mayorl 
(Il99ll) and Janson et al. (2012) we also used the probability 
density function for the orbital period and for the eccentricities 
(see Table O for G dwarfs, and the probability density function 
for the semi-major axis for M-dwarfs (see Table |4). The pro- 
jected distance R and the ARV were drawn performing a Monte 
Carlo simulation of Kepler's laws and its projection, consider- 
ing a uniformly random mean anomaly (M e [0, 27:]), inclina- 
tion (cos/ G [-1, 1]), longitude of the node (Q e [0,2;:]), and 
argument of the periapsis (oj e [0, 27i]). 

To compute the luminosity of the stellar companion given 
the mass of the star, we used relation l2l from iDurici (i2QQ3i) : 



x4.0 



t^it)'^ ifM>0.43Me; 
^-0-23(£r, ifM<0.43Me. 



(2) 



5.4.2. Statistical analysis of fortuitous alignments 

For the fortuitous alignment we are only interested in contami- 
nant stars whose flux can enter the fiber. With a seeing of 0.93, 
which is the mean value for the La Silla observatory, we mea- 
sured that the fraction of flux that enters the fiber of a star at 2.6 
arcsec from the fiber center is 2 x 10"^. If target and contaminant 
star have the same magnitude, this corresponds to an eff'ective 
Am of 19 within the fiber, and so its impact on the target RV can 
be considered null for our purpose. We then created a probability 
density function to have a contaminant star of a certain spectral 
type and magnitude within a radius of up to 2.6 arcsec around 
the target star, using the density tables from Appendix lAl 

For background/foreground stars within the 2.6 arcsec ra- 
dius we assumed a uniform distribution, and if a contaminant 



is present, the probability that the star is standing at a distance r 
from the target is proportional to r^, and the probability to have 
the fortuitous contaminant star between r and r -\- dr is given by 
equation [3l 



P(r, r + dr) = nilrdr + dr). 



(3) 



in which the members are self-explanatory. 

To calculate ARV we used the probability density function 
for the systemic RV, which was computed by fitting a Gaussian 
to the systemic velocity of stars from HARPS. This yielded an 
average and a dispersion of (yu, cr) = (7.87, 29.87). 

5.4.3. Statistical analysis results 

The probability to have a given impact on the RV calculation for 
each case is presented in Fig.|9l We can see that the most proba- 
ble scenario is to have a contamination of less than 10 cm/s. But 
contamination may happen and is more probable for far-away 
target stars, either because they become fainter, or because, in 
the case of a binary star, it is more probable that the compan- 
ion star flux enters the fiber. We can examine the target star 
G8, for instance. In this case, if the target star is 100 pc away, 
there is only a probability of 67% for a contamination below 10 
cm/s, and 14% probability for a contamination between 10 and 
100 m/s. There is a peak around a contamination of 10^ cm/s, 
which is, very probably, caused by physical binary contamina- 
tion. This level of contamination needs to be taken into account 
when choosing targets for planet searches. This becomes even 
more important when choosing target stars to be observed with 
ESPRESSO. 



6. Discussion 

We showed that we should not neglect the possibility that a stel- 
lar companion contaminates the RV calculations. This is spe- 
cially true for far-away stars, which have a higher probability 
to have two stars within the fiber. Nevertheless, we stress some 
assumptions we made that may have some impact on our work. 
Because spectra with high SN are scarce for late-M stars, in our 
statistical analysis we used the Ml spectra to represent the [MO, 
M5[ bin and the M3 spectra for the [M5, M9]. We also assumed 
that the star will have at most one stellar companion, because 
we found that the probability to have two contaminant stars is 
too low. We also did not take into account the eff'ect of the sky 
brightness if there was no stellar companion. 

We also stress that the Besan^on model is a probabilistic 
model, and thus is not free of errors. We consider it reasonable 
that the distribution of fortuitous alignments is a Poisson distri- 
bution of stars. Thus, the associated error is equal to Va^, where 
N is the number of stars per area element. 

The results presented in Fig. [7] may lead one to think that 
contamination may be stronger than it actually is, but we should 
not forget that this figure only shows the maximum impact 
of contamination as a function of Am, which occurs when 
\RV{contaminant)\ ~ 0.6 x FWHM(target). We also consid- 
ered that both stars are well-centered in the fiber and that all 
flux enters the fiber. This justifies our statistical analysis for the 
impact of the stellar companion within the fiber. With our statis- 
tical analysis we are able to conclude that it is important to take 
into account the possibility of contamination by a stellar com- 
panion within the fiber, especially when the RV precision enters 
the domain of the cm/s. Although the probability to have a con- 
tamination lower than 10 cm/s is above 50% for every target star 
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considered, we should pay special attention to far-away target 
stars, for which the impact of contamination can be higher than 
10cm/ s for 30 - 40% of the cases. 

6.1. Impact on typical cases 

As a consequence of the results presented in Fig IT] we present 
in Table [8] a study on the limit Am for which the maximum im- 
pact on the RV calculation is higher than 10 cm/s, 1, 10, and 
100 m/s. From this table we can see that for a G star with a 
contaminant M star, an impact higher than 100 m/s only occurs 
if Am < 2. However, in Table [5] from Sect. 12.21 we have seen 
that the probability to have a contaminant M star of my < 17 is 
very low. Thus, if we are only interested in target stars brighter 
than mv = 15, an impact higher than 100 m/s is very unlikely. 
We can also see that when observing an M star, an F star can 
never cause an impact higher than 100 m/s, and an impact of 
tens m/s is only possible if target and contaminant have almost 
the same magnitude (Am < 1.5). On the other hand, even for 
actual precisions, we may reach a potentially detectable contam- 
ination due to a Am = 9. In many cases the induced effect is not 
detectable because it is a constant effect. Nevertheless, there are 
conditions (variable seeing, centering and/or pointing problems, 
etc.) in which the induced effect will not be constant. Even so, 
a very particular configuration is necessary (and is highly un- 
likely) for the effect to create a periodic signal. A contamination 
stronger than 1 m/s on the RV calculation of a KO target star 
is possible if Am < 9. If the KO star is at 10 pc, it will have a 
my ~ 6, which means that contaminant stars with 6 <v< 15 can 
cause such a contamination. The probability to have one con- 
taminant star with 6 < my < 15 within the 0.5" radius fiber is 
0.02% only. 

These tabled effects correspond to the maximum case, i.e., 
to the difference between contamination and no contamination. 
Discussed cases are extreme in the sense that they are contam- 
ination and no contamination. Nevertheless, to extensively ana- 
lyze a variable contamination impact we would need to perform 
many hypotheses with uncertain parameters, such as seeing or 
flux counts. 

6.1 .1 . Sky brightness 

In our analysis of the sky brightness impact we assumed that sky 
brightness is mainly caused by the moon, and so its spectra will 
be that of a star with spectral type G2. The new-moon sky bright- 
ness should then be considered as a lower limit for the contami- 
nation from the sky brightness with a G2 spectrum. However, we 
should not forget that sky brightness has other sources and may 
depend on the sky zone. Sky spectra from the Galactic disk are 
surely different from the Galactic halo. Moreover, there is also 
the possibility to have a galaxy as contaminant, which was not 
considered. Nevertheless, we consider it a good approximation, 
or at least a good starting point, to treat the sky brightness as a 
G2 star. 

If observations are carried out on a night with moon, and if 
cirrus are present, the contamination by the moon will be en- 
hanced. 

iTripathi et al.l (l2010l) . reported two outlier RV data points 
caused by cirrus during their observations of WASP-3, an F7 
star of my = 10.7, on a full-moon night. Their RV measurements 
present a redshift of 140 m/s and of 49 m/s with respect to their 
best-fitting models. If we assume that cirrus will reflect some of 
the moonlight, we have a G2 contaminant. The maximum value 



Table 8: Limit Am for which the maximum impact (Imax) on the 
target RV is higher than 0.1, 1, 10, and 100 m/s, for each combi- 
nation of target spctral type (TST) and contaminant spectral type 
(GST). 



CST 

TST W[m/s] F2 F5 G2 G8 KO K5 Ml M3 
G2 



G8 



KO 



K5 



Ml 



M3 



considered for all practical effects for the visual magnitude of 
this contaminant will be that of the target star (if it were higher 
than the target magnitude, the target star would be no longer 
visible). Thus, for a G2 target star with my = 10.7 and a contam- 
inant with the same spectral type and magnitude, the maximum 
contamination will be 2.7 Km/s. This corresponds to 0.15% of 
the moonlight flux, considering that 100% of the moonlight flux 
corresponds to a star of magnitude 3.65. If 0.01% of the moon- 
light reaches the fiber, corresponding to a contaminant of spec- 
tral type G2 with my = 13.7, the contamination would be 140 
m/s. For a contamination of 49 m/s, 0.003% of the moonlight 
flux would be enough, co rresponding to a G2 contaminant with 
my = 14.8. Because the ITripathi et al.l (l2010l) data were taken 
with the High Resolution Echelle Spectrometer (HIRES) on the 
Keck I telescope, and we used HARPS spectra in our simula- 
tion, and also because we did not calculate the impact for a tar- 
get of spectral type F, these contamination values are presented 
here just as a reference, since we cannot directly compare them. 
Nevertheless, we consider that their hypothesis of cirrus contam- 
ination is valid. 

The sky brightness contamination is not easily removed. 
Although its signal may be strong enough to contaminate the RV 
calculations, it will be too weak to be adequately characterized, 
fitted, and removed, even for a long- slit spectrograph. 

6.1.2. Kepler stars 

Our work can also be a valuable asset for the follow-up of pro- 
grams such as Kep ler. Kepler stars have 11 < my < 18 (see 
iBatalha et al.l [2012h . and consequently, measuring the RV be- 
comes challenging. Beyond the problem of these being faint 
stars, there is also the problem of the false-positive planet 
transit. As an example, we consider the case of Kepler- 14b 
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(iBuchhave et al.|[20Tn) . a planet orbiting an F star with a com- 
panion star of nearly the same magnitude and only 0.3 arcesec 
of sky-projected angular separation. Buchhave and collaborators 
were not able to detect a second peak on the CCF and concluded 
that the two stars probably have nearly the same RV. They com- 
puted an impact of 280 m/s on the RV caused by the companion 
star. A direct comparison with their work is not possible, because 
their observations were carried out with the Fiber-fed Echelle 
Spectrograph (FIES) at the 2.5m NOT at La Palma, and with the 
High Resolution Echelle Spectrometer (HIRES) mounted on the 
Keck I on Mauna Kea, Hawaii. Moreover, our simulation only 
contemplates target stars of spectral GKM, and this is an F star. 
Still, performing a rough comparison, we see that our simula- 
tions indicate an impact of ~ 276 m/s for a G2 target star with a 
contaminant of the same spectral type. Am = 1 , and ARV = 1 . 

7. Conclusions 

This study on the impact of stellar companions within the fiber 
of RV calculations allows us to conclude that we should not ne- 
glect the possibility of contaminant flux from stellar companions 
stars, specially if we are observing far-away stars. On average, 
if we are observing a G or a K star, the contamination may be 
higher than 10 cm/s if the difl'erence between target and contam- 
inant magnitude is Am < 10, and higher than 1 m/s if Am < 8. 
If the target star is a star of spectral type M, a Am < 8 is enough 
to obtain the same contamination of 10 cm/s, and a contami- 
nation may be higher than 1 m/s if Am < 6. We also showed 
that sky brightness should not be discarded, particularly on full- 
moon nights, if one observes faint target stars (m^ > 10). 

These results will allow us to more wisely choose tar- 
get stars to be observed with instruments such as ESPRESSO 
and CODEX, and they provide reference values for the difl'er- 
ent cases of contamination possible in fiber-fed high-resolution 
spectrographs. We also stress that there are diagnosis methods 
that should be capable of detecting most of the blends that 
mimick planets. We can detect a blend through bisector anal- 
ysis of the C CF or from correlations using diff'erent templates 
(ISantos et al.ll2002h . 
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Fig. 7: Maximum impact of \Am\ on the \RVCt\ for a target star of spectral type a) G2, b) G8, c) KO, d) K5, e) Ml, and f) M3, and 
contaminants of spectral types F2 (star), F5 (hexagon), G2 (squares), G8 (diamonds), Kl (circles), K5 (triangles). Ml ( x crosses), 
and M3 (+ crosses). 
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Fig. 9: Close-up of the distribution of the expected contamination , \RVCt - RVC^l, on the radial velocity, |7?yQ|, of a target star of 
spectral type a) G2, b) G8, c) KO, d) K5, e) Ml, and f) M3, and at distance d, caused by a baclcground or a companion star. The 
green thiclc solid line shows the distribution if the target star stands at a distance of d = 1 pc, the blue dotted line for d = 5 pc, blaclc 
dashed line for d = 10 pc, magenta dash-dotted line for d = 50 pc, and red solid line for 100 pc. Numbers in each bin correspond to 
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Fig. A.l: Number of stars in each bin of spectral type and magnitude for an area of the sky with a radius of 3°. 
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Fig. A.2: Density of stars in each bin of spectral type and magnitude for an area of the sky with a radius of 0.5 arcsec. 
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Appendix B: Method details 

To control the magnitude of the contaminant star we first used the relation IBJI to calculate in the target reference frame, 

mt-mc = -2.5 x log(F,/F,) , (B.l) 

where Ft and Fc are the fluxes of the target and of the contaminant. The magnitudes are now calibrated and can be changed by 
multiplying the flux by a factor of 

/ = 10-^/2-^ (B.2) 

where x is the difl'erence between the new chosen magnitude and the true magnitude of the star in the target-star system. 

To change the RV of the contaminant spectra we also had to modify them. The first challenge when processing our co-added 
spectra with the HARPS RV pipeline is that it will take into account the header (and proprieties) of the parent spectra. To compensate 
for this we first need to replace the header of the contaminant spectrum file with that from the target, and then run the cross correlation 
function (CCF) recipe of the RV pipeline. The next step in this process was to shift the contaminant spectra. To assess by how much 
the contaminant spectra had to be shifted, we used equation lB3l 

^ I iRVCc - RVCt - ARV\\ 
Ac = At xll + l , (B.3) 



where Ac are the new wavelengths to derive the desired ARV (difl'erence between the RV of the target and the contaminant). At 
are the wavelengths from the wave file corresponding to the target, and c is the light speed in vacuum. These new wavelengths 
were used to calculate, by linear interpolation, the number of pixels by which the contaminant needed to be shifted. This gives a 
new XX axis that it was used, recurring to another linear interpolation, to shift the contaminant spectra. With the target RV for the 
contaminant, the next step is the sum of the target with the modified contaminant spectra, which is the sum of their fluxes. The 
shift of the contaminant spectra should not be larger that the maximum barycentric Earth radial velocity (BERVMX) of the star, 
otherwise there is the risk that the spectral lines enter and leave the orders when performing the CCF calculation. BERVMX defines 
the excluded spectral lines from the correlation due to the orbital movement of the Earth toward the star. Our goal is to introduce a 
shift equivalent to a BERV value lower than the |BERVMX|, so that we do not need to exclude new lines. 
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